Verwey transition in Fe 3 04 at high pressure: quantum critical behavior at the onset of 

metallization 
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We provide evidence for the existence of a quantum critical point at the metallization of magnetite 
Fe304 at an applied pressure of p c ~ 8 GPa. We show that the present ac magnetic susceptibility 
data support earlier resistivity data. The Verwey temperature scales with pressure TV ~ (1— p/pc)", 
with v ~ 1/3. The resistivity data shows a temperature dependence p(T) — po + AT", with n ~ 3 
above and 2.5 at the critical pressure, respectively. This difference in n with pressure is a sign of 
critical behavior at p c . The magnetic susceptibility is smooth near the critical pressure, both at 
the Verwey transition and near the ferroelectric anomaly. A comparison with the critical behavior 
observed in the Mott-Hubbard and related systems is made. 
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Magnetite (Fe304) is among the oldest and most fasci- 
nating materials throughout the history of science. It is 
ferrimagnetic below the Curie temperature Tq ~ 860 K 
and semiconducting with a room temperature resistiv- 
ity of p(T = 300 K) ~ 10- 2 ficm. The stoichiometric 
material undergoes a discontinuous Verwey transition 1 
to a higher-resistivity phase at a temperature of Ty ~ 
(121 ± 1) K, at which the resistivity at ambient pressure 
jumps up by about two orders of magnitude upon cool- 
ing. The transition is of a semiconductor-semiconductor 
type at lower pressures, although the exact difference be- 
tween them remains somewhat elusive. In addition, the 
net magnetic moment has the value p ~ 4.1/^s per for- 
mula and approximately corresponds to the spin 5 = 2 
in Fe 2+ high-spin configuration in the octahedral position 
(the remaining two Fe 3+ ions, one in tetrahedral and one 
in the octahedral positions with spin 5 — 5/2 have their 
moments antialigned). The high-spin state is induced by 
Hund's rule coupling among 3d electrons. 

Essential progress was made when samples with an 
extremely well controlled composition were synthesized 2 
and studied.— Notably, the change of the Verwey transi- 
tion from first- to second-order has been discovered, de- 
pending on the off-stoichiometry parameter S of Fes-aO^. 
Also, the metallization of Fe3C>4 has been achieved^ at a 
pressure of p c ~ 8 GPa. 

The primary purpose of this work is to systematize the 
change in the properties of Fe304 with pressure and to 
single out universal features of these results. During the 
course of this work we reanalyzed the earlier resistivity 
data* and have discovered that the combined magnetic 
susceptibility and resistivity data point to the possibil- 
ity of a quantum critical point existing at the metalliza- 
tion threshold. This is explicitly verified by determining 
the pressure dependence of the critical exponents of both 
the transition temperature Ty (p) and the resistivity just 
above the transition p(Ty + ), as well as that of p(T) at 



the critical pressure and above it. Also, the pressure de- 
pendence of the two low-temperature anomalies in the 
ac susceptibility, presumably connected with the appear- 
ance of ferroelectricity, evolve with increasing pressure, 
but this aspect of the work will not be detailed here. 

The importance of the presence of a quantum-critical- 
point (QCP) at the localization- derealization bound- 
ary cannot be underestimated as there are no well de- 
fined universality classes as yet for this electronic tran- 
sition. Previously, a classical critical point (of the Van 
der Waals-type) was determined^ in a three-dimensional 
system Cr-doped V2O3, at T cr u — 457.5 K and p cr it — 
3.7 kbar, the existence of which was suggested earlier £ 
The critical behavior in this case has the critical expo- 
nents of the mean-field type, as the T cr n value is rather 
high. Recently, a critical point has been seen 7 in a 
quasi-two-dimensional organic solid at T cr n — 39.7 K 
and p c — 24.8 MPa, with completely different critical 
exponents than those for the above three-dimensional 
case. A true quantum critical point (at T = 0) at the 
localization-delocalization boundary has been detected 
at the Kondo semiconductor - non-Fermi liquid bound- 
ary upon chemical substitution in CeRhSbi-^Sn^i^ In 
the last system, it is connected with the delocalization of 
4/ electrons due to Ce 3+ and is induced by the change 
in carrier concentration upon doping. We suggest, that 
the present system (Fe304) represents the first exam- 
ple of a system with a true quantum critical point (with 
T C rit = 0) in a monocrystalline and stoichiometric sys- 
tem in a spin-polarized state, both in the semiconducting 
and the metallic states. In contrast, a different type of 
behavior has been observed^S in the antiferromagnetic 
system NiS2- a; Se :r , in which the Hund's rule coupling, 
leading to the high-spin 5 = 1 state in this case, also 
occurs. The differences are caused by a disparate mag- 
netic ordering, as well as by the circumstance that, as we 
show, the metallic state of Fe304 is that of a non-Fermi 
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FIG. 1: Pressure dependence of ac magnetic susceptibility (at 
a frequency of 120Hz) for the stoichiometric single crystal; 
note different T scales for the real (x') and the imaginary 
(x") parts of x- Right inset: illustration of the definition of 
Tv from maximum of dx' /dT when the sharp step in x' was 
not observed (i.e. forp > 3.5 GPa). Left inset: T max indicates 
the position of the anomaly near 50 K, whereas T m i n the lower 
anomaly (indicated by the vertical arrow). 



liquid, whereas that of NiS2-xSe x is an almost localized 
antiferromagnetic Fermi liquid. 

A single crystalline magnetic sample used for the ac 
magnetic susceptibility measurements, was cut from a 
larger crystal grown from the melt using the cold crucible 
technique (skull melter)^, at Purdue University. The 
crystal was then annealed under CO/CO2 gas mixtures 
to establish the stoichiometry-2 and rapidly quenched to 
room temperature to freeze in the high temperature ther- 
modynamic equilibrium. Although this procedure gen- 
erates octahedral defects, most of the low temperature 
electronic processes remain intact, as is evidenced by the 
sharp Verwey transition temperature Ty and the exis- 
tence of the low tempearature anomaly (observed only in 
the best crystals). 

In Fig. 1 we display the temperature dependence of 
the ac susceptibility X = x' + *x" components for se- 
lected pressures. The upper panel provides the real part 
x', for which three characteristic temperatures are: The 
Verwey temperature Ty (the uppermost jump), and the 
two temperatures associated with ferroelectricityr^ for 
Tmax ~ 60 — 50 K, and the lower anomaly at T m i n ~ 
20 K. With increasing pressure, all of them shift towards 
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FIG. 2: Pressure dependence of the Verwey transition tem- 
perature Tv- We mark the data points from this work as 
solid squares (blue online) and those of Mori et al.— by open 
squares. Solid line represents the fit specified and ending at 
quantum critical point (QCP). Inset: same in the doubly log- 
arithmic scale. Open triangles show the respective critical 
behavior for slightly nonstochiometric (S = 0.007) sample. M 
denotes a stable metallic state, whereas SC labels the semi- 
conducting state below Ty. 



lower temperatures and eventually disappear for pres- 
sures p > 7.5 GPa. Note that pronounced losses are 
observed near T max . They are most probably due to the 
ferroelectric transition and associated with it ferroelectric 
domain dynamics. In the insets we mark the tempera- 
tures Ty, T max , and T m i n near the critical pressure. The 
first two temperatures approach each other, as well as 
decrease with increasing pressure. 

To explicitly demonstrate the evolution of the Verwey 
transition with pressure, we plot (Fig. 2) the Ty(p) de- 
pendence by combining our magnetic susceptibility re- 
sults and the resistivity data taken from Ref. |j. In the 
inset we redraw the data on a doubly logarithmic scale 
for both the almost perfectly stoichiometric (Fe304.oo2) 
and the slightly nonstoichiometric (Fe304.oo7) samples 
to demonstrate the exponential fit quality, as well as to 
show the difference between the two cases. We observe 
a clear exponential behavior of Ty(p) = A{\ — p/p c ) v 
above 3.5 GPa for the stoichiometric sample, with the 
critical pressure p c ~ 8 GPa and the value of the ex- 
ponent v ~ 1/3. Actually, the p(T) data of Mori et al4 
already shows metallic conductivity at p = 7.9 GPa down 
to about 2 K, but a percent inaccuracy in determining 
p c is acceptable given the inaccuracy of determining Ty 
for p > 7 GPa. As the data in the inset shows, both p c 
and v depend strongly on the sample stoichiometry. For 
example, for S = 0.007 we have p c ~ 6.1 GPa and the 
critical exponent is reduced to the value v ~ 1/6. This 
means that the nonstoichiometry is not equivalent to a 
reduction of critical pressure. 

The nature of the metallic state (M) above p c should 
be characterized in greater detail. This is particularly 
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FIG. 3: Pressure dependence of the electrical resistivity just 
above the Verwey transition 4 (at T = Tv+). The residual 
resistivity obtained from the fit is po — within the fitting 
error. 
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FIG. 4: Temperature dependence of electrical resistivity 
p(T) = po + AT n of Fe 3 O 4 .002 at p c = 8 GPa. The solid line 
(red online) corresponds to po = 35.5 pQcm and n = 2.52. 
Open symbols show the sensitivity of the fits to po value: 
slight variations of the po value do not appreciably change 
the fitting accuracy. 



important because for p < 7.5 GPa the state above Ty 
is clearly a semiconductor. A similar type of behavior to 
that exhibited in Fig. 2, though without any reference to 
quantum critical behavior, was suggested by Rosenberg 
et alJ^ Also, the magnetic susceptibility and the resistiv- 
ity discontinuities in Fig. 2 follow the same dependence 
Tv(p)- 

To provide quantitative evidence that we observe a 
criticality near the T cr n — 0, p — p c point, one re- 
quires a proper definition of the order parameter, as well 
as its dependence on p — p c . In the previous studies 
of criticality near the insulator-metal transition^ the 
role of the order parameter was played by the conduc- 
tivity difference a(p) — cr(p c ) or equivalently, directly 
by the resistivity^ The two transport quantities mimic 



the essential difference between carrier concentrations on 
both sides of the electronic transition in the equilibrium 
state. Additional evidence is provided by the behavior 
of physical quantities near p c . For that purpose, our 
data from Fig. I is not directly useful, as the suscepti- 
bility is continuous upon approaching the system metal- 
lization. This in turn means, that magnetism is robust 
under the applied pressure, as no clear signature of crit- 
ical magnetic behavior has been observed up to 300 K. 
In view of that circumstance, we reanalyzed the data 
of Mori et al4^ and plot in Figs. 3 and 4, respectively, 
the resistivity p(Ty+) just above the Verwey transition 
as a function of p — p c , as well as p(T) at p c . We ob- 
serve that, in the former situation we have roughly the 
scaling p(p) — R[(p c — p)/p c ] a , with R ~ 117cm and 
a ~ 5/2, whereas at p = p c : p(T) = po + AT n , where 
po — 35.5/zficm, and n ~ 5/2. These results differ from 
those for Mott-Hubbard systems«i£ obtained when the 
critical temperature T cr i t > 0. However, it should be un- 
derlined again that, as only in the present case T CT n = 0, 
we are dealing with a true quantum critical regime for 

P^Pc- 

In Fig. 5 we plot the temperature dependence of the 
resistivity in the metallic phase (for p — 9 GPa). The 
resistivity p(T) is roughly proportional to T 3 and the 
value of the residual resistivity is considerably reduced. 
The reduction of po upon further metallization means 
that scattering on charged impurities (Fe 2+ or Fe 3+ va- 
cancies) is reduced, since they are screened more effec- 
tively with increasing pressure. Note also the extreme 
sensitivity of the fitting upon the choice of the po value. 
Parenthetically, this sensitivity may be regarded as an 
additional criterion of the sample quality. 

The well denned exponential behavior of both the 
pressure and temperature dependencies of the resis- 
tivity match neither those observed near the quan- 
tum critical point for a ferro- and antiferro-magnetic 
system,™ nor those mentioned here for the antiferro- 
magnetic systems near the critical point for the Mott- 
Hubbard transitions^^ In other words, the metalliza- 
tion of FC3O4 represents a separate universality class of 
phase transitions. In connection, it would be very in- 
teresting to compare them in detail with those^ for the 
spin-canted antiferromagnet MS2 , a Mott-Hubbard insu- 
lator at p = 0, for which the Hund's rule also plays as an 
important role as for the ferrimagnetic state of Fe3 04. 

A basic question of relating the observed metalliza- 
tion to the microscopic properties of the system re- 
mains. On phenomenological grounds, one can relate 
the thermodynamics of the transition to the charge- 
order transition in the octahedral sites, connected with 
the freezing of the resonating electron between Fe 3+ - 
Fe 2+ configurations of the two ions at Ty*^ Consider- 
ing the extra electron (Fe 2+ =Fe 3+ + e~) as undergoing 
Mott-Wigner localization,— one can write down a phe- 
nomenological Ginzburg-Landau functional, albeit with a 
fermionic form of the configurational entropy for the spin- 
less (i.e. strongly polarized) fermions, and rationalize a 
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FIG. 5: Temperature dependence of resistivity p(T) = po + 
AT n in the metallic (M) regime for Fe3C>4.oo2- The solid line 
corresponds to po = 7.25 pftem and n — 3.05. The fitting 
accuracy does not change under a minute change of po, as 
shown. 



changeover from the first-to second-order upon varying 
the system stoichiometryj^ 

In recent years the above picture was investigated using 
spectroscopic methods. 16 A pressure induced transforma- 
tion from the inverse to the normal spinel structure was 
claimed to occur XL This seems to be difficult to reconcile 
with the fact, that a transformation of this type should 
also take place at p = 0, which is not observed. In other 
words, the point Ty(p = 0) should become also a ter- 
minal point of the first-order line Ty(p). Nevertheless, 
in spite of the complications of the situation, in which 
both strong electronic correlations and electron-lattice 
coupling are important, a partial charge ordering was 
detected by X-ray scattering^ This last result means, 
that strong correlations may play a primary role in the 
transition. This is because a relatively simple exponen- 
tial behavior of Ty(p), p(T,p c ), and p(T,p > p c )) means, 
that the underlying mechanism of the transition should 
be relatively simple physically. Otherwise, the change 
of Ty from above 120 K to zero should be sizably in- 
fluenced by e.g. the difference in thermal occupation of 
the low-energy phonons and hence lead to nonuniversal 



behavior. 

Having in mind the underlying simplicity of the elec- 
tronic correlations as the primary force, we can easily 
rationalize the Ty ~ (p c — p) 1 ^ 3 dependence. Namely, 
regarding strongly correlated polarized 3d electrons as 
spinless fermions, we can write down the Landau free 
energy of the electronic system in the form^^ 

F = F + - k B T[^> ln(*/2) + (1 - *) ln(l - *)] 



Fq + a(T - T v )y + b^ 2 



(1) 



where E(^f) = £q + ei^t + t2^ 2 + ■ ■ ■ is the system en- 
ergy, the third term the relative fermionic entropy, \E' 
is the density of itinerant electrons (resonating Fc 2+ «-» 
Fe 3+ electrons in the octahedral sites), and e , . . . , b are 
constants. The coupling of the electron subsystem to 
the lattice introduces an additional term of the Holstein 
type AtfAx/a = X^{AV/V) 1/3 , where Ax/a and AV/V 
are respectively, the relative lattice-parameter and vol- 
ume changes under pressure, and A is the local electron- 
lattice coupling constant. Assuming that, ei ~ A(. . .) 
and since-^ AV^ ~ p c — p, the change of Ty induced by 
applied pressure is Ty = A'(p c — p) 1 ^ 3 , which is in agree- 
ment with what we observed. 

In summary, we have demonstrated, that the metalliza- 
tion of magnetite under pressure can be described by a 
simple exponential scaling of physical properties, includ- 
ing the transition temperature Ty, the resistivity near 
the critical pressure p(p, Ty + ) at the transition, and as a 
function of temperature, p(p = p c ,T), at p c and above. 
These properties provide direct evidence of critical be- 
havior with a quantum critical point at p c rj 8 GPa. It 
is also suggested, that the strong electronic correlations 
lead to a partial charge order (and concomitant with it 
partial localization) and may be the source of the ob- 
served volume change, as is the case for Mott-Hubbard 
systems. 
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